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A  New  Multiplier-Free  Structure  for  DFT  Based  on  Delta  Modulation 


Tang  Kaijun,  Siang  Jingcheng,  Huang  Shunji 
(Chengdu  Institute  of  Radio  Engineering) 

Abstract  A  new  implementation  of  the  DFT  based  on  delta 

modulation  (DM)  is  proposed.  It  is  different  from  the 
conventional  implementation  based  on  PCM.  The  hardware  structure 
of  the  method  is  constructed  without  A/D  converter  and 
multiplier,  so  it  is  simpler  and  of  lower  cost.  Its  SNR  can 
approach  that  of  the  conventional  one,  such  as  FFT.  Computer 
simulations  demonstrate  that  for  deterministic  signals,  results 
obtained  agree  well  with  theoretical  analyses,  for  band-limited 
Gaussian  signals;  if  step  size  is  suitably  selected,  results 
similar  to  the  conventional  DFT  can  still  be  obtained.  Because  of 
the  limitation  of  the  working  speed  of  the  device,  this  method 
will  be  suitable  for  processing  speech  and  seismic  signals. 
Finally,  the  hardware  structure  which  simply  consists  of  a  ROM, 
an  adder  and  some  other  auxiliary  circuits  is  given. 

(Manuscript  received  on  December  30,  1984  and  revised  version 
received  on  June  20,  1986.) 

(Research  sponsored  by  National  Science  Foundation  of  Chinese 
Science  Academy.) 

I.  Introduction 
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Delta  modulation  (DM)  is  a  simple  coding  method  for  various 
signals.  In  the  past  DM  was  mainly  used  in  coding  signals  for 
transmission.  Recently,  application  of  DM  has  been  extended  to 
processing  of  digital  signal  and  a  family  of  DM  wave  filters [2-7] 
have  been  proposed.  However,  research  aimed  at  applying  DM  to 
spectrum  analysis  is  rare.  The  purpose  of  this  paper  i3  to 
utilize  the  relationship  between  DM  serial  spectrum  and  PCM 
serial  spectrum  to  carry  out  spectrum  analysis  by  DM  method. 

II.  Delta  Modulation  Series  Fourier  Spectrum 

Binary  DM  series  is  considered  in  this  paper.  Assuming  that 
proper  step  size  and  sampling  frequency  were  selected,  the 
probability  for  the  slope  to  be  overload  is  very  small  and  can  be 
ignored. 

A  DM  series  is  defined  as: 

dx(n')  ^  sign(x(n)  —  &(n  —  I)),  (1) 

A 

where  x(n-l)  is  the  predicted  value  at  n-1  instant.  Because  only 
linear  DM  system  is  studied  in  this  work,  the  step  size  A  is  a 
constant.  Figure  1  shows  the  delta  modulator  and  figure  2  is  the 
equivalent  liberalized  block  diagram  of  figure  1. 

Based  on  figure  2  we  have: 

ddx(ti)  —  x(f»)  —  *(»  -  1)  +  —  1),  (2) 

where  eq(n)  is  the  quantified  error  series.  Statistical  analysis 
shows  that  it  can  be  considered  as  the  white  noise  homogeneously 
distributed  between  [-  ^  ,  A ]  . 

The  DM  series  Fourier  spectrum  can  be  defined  by: 
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Fig.  1:  Delta  modulator 


Fig.  2:  Linear  equivalent  diagram  of  delta  modulator 
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N-l 

4  —  0 ,  1 , •  •  • ,  .v  —  1 ;  (3) 

•-0 

or 

.V-  1 

*Ak)  “  2  dx(n)h{N  —  1  —  n,  *),  {  —  0,  1,- •  • ,  N  —  1;  (4) 

where  h(n,k)-Wj|~^^n+^)  . 

Figure  3  shows  the  circuit  diagram  of  analyzed  DM  spectrum 
and  figure  4  is  the  linear  equivalent  diagram  of  the  entire 
system.  Equation  (4)  is  actually  the  roll  product  of  the  input 
series  dx(n)  and  the  pulse  response  function  h(n,k)  of  FIR  wave 
filter. 

The  relationship  between  the  Fourier  spectrum  Xq(k)  of  pulse 
coding  modulation  series  and  the  Fourier  spectrum  Xj(k)  of  DM 
series  is: 


Xt.U)  -  *«.-,<*)  +  ax  MU  (5) 

where  the  definition  of  Xqn(k)  and  Xjn(k)  are  defined  by: 

N-l 

+  *-0,  -l;  (6  ) 

rm  0 

-  N-| 

&  2  dx(-n  +  r),yk>  4  -  0.  *»•••,  N  —  1.  (7) 

'  -  0 

and  xq(n)  in  equation  (6)  is  the  quantitative  magnitude  of  x(n) . 

Let  xq(n)-x(n)  and  from  the  relationship  x(n)-x(n-l)+  4dx(n) 
we  can  obtain  equation  (5) .  Equation  (5)  shows  that  the  Fourier 
spectrum  of  PCM  can  be  obtained  by  summing  the  Fourier  spectra  of 
DM. 

The  sampling  rate  parameter  K  is  defined  as  T/^t,  where  T  is 
the  Niquest  sampling  spacing  and  At  is  the  DM  sampling  spacing 
and  N=KM  where  M  is  the  number  of  sampling  points  based  on  the 
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Fig.  4t:  Linear  equivalent  diagram  of  the  entire  system 


Niquest  sampling  criterion.  Because  the  signal  band  width  is 
fixed,  after  adopting  the  DM  system  the  sampling  rate  is  enhanced 
by  K  times  and  X(k)-0  within  one  period  of  the  cycle  domain  when 
M/2<-k<N-(M/2) . 

After  adopting  DM,  the  speed  for  the  corresponding  addition 
operation  is  improved  by  K  times  and  the  multiplication  operation 
of  the  entire  system  are  eliminated.  For  signals  with  narrow  band 
width  such  as  speech  signals,  the  speed  for  addition  operation 
will  provide  a  healthy  margin  and  this  plan  will  work  fine. 


III.  Quantitative  Error  Analysis  of  Delta  Modulation 

The  transformation  from  the  non-linear  DM  system  to  linear 
system  is  by  the  introduction  of  quantitative  noise.  It  has  been 
discussed  above  that  the  treatment  of  quantified  noise  as 
homogeneously  distributed  white  noise  is  reasonable.  This 
treatment  also  simplifies  the  system  analysis. 

Assuming  the  input  is  the  white  noise  with  power  cr^2  then 
the  power  of  the  output  signal,  based  on  figure  4,  is-. 

£[|X(*)n-Afci,  (8) 

Similarly,  the  power  of  the  output  noise  is: 

(9) 

Therefore,  the  output  signal  to  noise  ratio  is: 


£[|X(t)H  _  gj 
Et|£(*)l,j  < 


(10) 


where  Giloi,  —  CK'*,  is  the  relationship  between  the  input  signal 
to  noise  ratio  and  the  sampling  rate  coefficient  K  under  the 
condition  where  there  is  no  slope  overload  and  maximum  signal  to 
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noise  ratio, (7]  C  is  a  proportionality  constant.  If  the  input  is 
the  band-limited  Gaussian  signal,  then  C-6/(aTr)2  and  a  should  be 
determined  by  the  probability  of  occurrence  of  slope  overload. [7] 
In  PCM  systems,  the  increase  in  the  input  signal-to-noise 
ratio  is  accomplished  through  extending  the  word  length  of  A/D 
convertor,  while  in  DM  systems,  it  is  accomplished  through 
increasing  the  sampling  coefficient.  Therefore,  in  digital  signal 
processing,  DM  method  can  be  used  to  replace  PCM  methods  and 
achieve  the  same  result.  Hence,  DM  method  is  an  effective  method 
in  realizing  digital  signal  processing. 

(In  reference  (71,  -  ex*,  and  the  assumption  * - -£L.  was 

made  by  its  author.  If  the  noise  in  the  band  width  of  the  signal 
is  considered  we  would  obtain  and  a’jaitmCK>  .) 

IV.  Analysis  of  Limited  Word  Length  Effect 

There  are  two  realization  methods  in  digital  signal 
processing;  fixed  point  realization  and  floating  point 
realization.  Because  of  the  limitation  of  word  length  of 
available  hardware,  only  limited  number  of  digits  can  be  used  and 
certain  error  is  brought  to  the  actual  result.  The  following 
analysis  shows  the  effect  of  these  two  realization  methods  on 
system  function. 

1.  Fixed  point  realization  limited  word  length  effect 
As  shown  in  figure  4,  because  of  the  absence  of 
multiplication  operation  there  is  only  one  kind  of  error;  namely. 
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the  coefficient  quantization  error.  As  discussed  in 
1 iterature  14]  this  is  a  homogeneously  distributed  white  noise. 

However,  in  order  to  ensure  no  overflow  when  summation  is 
taken,  the  summation  was  carried  out  in  groups  as  shown  in  figure 
5. 

Based  on  equation  (5),  the  relationship  between  X(k)  (the 
subscripts  q  and  n  are  deleted  and  X(k)  is  treated  as  the  Fourier 
spectrum  without  coefficient  quantization  beginning  with  n=0)  and 
X^(k),  note  that  0<k<^/2  and  K>“8[9]  and  1 2ttk/N  l<rr/K<  1  if  the 
input  signal-to-noise  ratio  is  to  be  greater  than  20dB,  is 

X(t)~ 


To  ensure  the  absolute  magnitude  of  the  quantity  in  the 
bracket  in  equation  (11)  is  smaller  than  1,  it  is  multiplied  by  a 
constant  C  and  let 

*)  -  w**'*0- 

1 

with  a  suitably  chosen  constant  C,  hi(n.A)  can  be  sufficiently 
expressed  in  terms  of  words  with  limited  length  as  figure  5 
shows . 

Another  error  is  introduced  by  this  realization  method;  the 
digit-shifting  operation  error,  which  is  treated  again  as  a 
homogeneously  distributed  white  noise  and  can  be  expressed  by 
q(p,n)  (p-1,2,  ...,  m»log2N;  n-0, 1 , . . . ,N-1 ;  where  N  is  the  powers 
of  2).  Its  real  and  imaginary  parts  are  taken  as  -2“'b  D,  0, 
2~(b+l)  With  equal  probability.  From  figure  5, 

X({)  -  X  (i-Y°t,N  r(N  -  1  -  »)C*(»,  *),  ( 12) 
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V—  | 

*«(*)  -  X  (±)>Ct,N  x(N-  1  -  n)CA,(»,  *),  (13) 

where  Xq(k)  is  the  Fourier  spectrum  with  quantized  coefficients 

The  error  of  coefficient  quantization  is 
ECO  -X(0  -*«(*) 

X  r(N  -  1  -  »)Cr(»,  O 


■  *  0 

N’-  1 


-  X-*('v- 1 

icQ  iV 

The  power  of  the  output  noise  is: 

CVX 


E(!E(On 


iV 


(14) 


(15) 


where  crc2«2-2b/6K,  and  crx2  is  the  power  of  input  signal  . 
Similarly,  if  the  input  is  the  white  noise  signal,  the  output 
signal-to-noise  ratio  is 


El\E(k)\*]  ,  2^ 

E[|X(OI*]  '  6  K’ 


(16) 


Similar  to  the  above  analysis,  the  error  corresponding  to 
any  digit-shifting  of  the  s  level  is 


N/j'-'-l 


X  (~)  *  ~~  1  “  «) 


•,g  '  2 

2, 


—  X  —  *iN  —  1  —  n)v(f>  "). 

i.g  N 


(17) 


The  total  noise  power 


Et|E,(OI*]-EUff.<*.  l)|’]  +  E[|E.(*f  2)|»]  +  ...  +  E[|E,(*ti)|»] 


“  —  triad  +  —  critri  +  •  •  •  +  2tri<ri 
N  N 

4«rJ<ri. 

and  the  output  signal-to-noise  ratio  is 


E[|E,(«n  _  ,,,  ,  _  4 


(18) 


(19) 


where 


«ri-  2(-i-2-K*+0  +  -l-2-,<‘+0'\--i-  2-J».  © 

'3K  3K  /  3  K 
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Comparing  equations  (19)  and  (20),  it  is  obvious  that  the 
error  from  digit-shifting  is  greater  than  the  input  of 
coefficient.  Therefore,  digit-shifting  noise  is  the  major  source 
of  noise.  This  result  is  consistent  with  the  noise  signal  of 
fixed  point  realization  when  FFT  is  calculated. 

2.  Floating  point  realization  limited  word  length  effect 

Even  though  fixed  point  realization  is  easier,  it  requires 
longer  word  length  and  more  adders.  The  floating  point 
realization  method  discussed  below  will  avoid  this  problem. 

Figure  6  is  the  equivalent  statistical  block  diagram. 

Similar  to  the  above  analysis,  the  coefficient  quantization 
error  e(n)  and  operation  error  qi(n)  are  treated  as  homogeneously 
distributed  white  noise  with  a  average  value  of  0  and  square 
difference  of  2/3K2_2b<?  and  2/3j^2_2br,  where  bq  and  br  are  the 
coefficient  quantization  word  length  and  operation  round-off  word 
length,  respectively. 

Similar  to  the  analysis  in  reference  [1],  if  the  input  were 
the  white  noise  signal,  then  the  output  signal-to-noise  ratio  is 

E[|E(OI*]  ^  _2_  2~Ur 

E[|*(t)n  IK  3  '  (20) 

Computer  simulation  has  shown  that,  if  the  input  signal-to- 
noise  ratio  is  40dB  then  K  has  the  value  of  32.(91  If  M  is  taken 
as  128,  then  N  is  4096.  Because  generally  br  >>  bq  and  if  the 
latter  is  ignored,  then  the  above  equation  can  be  approximated  by 

i 

E[|E«)n  ^  2  ,,f  (21) 

E[|*U)I’J  3K 
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Fig 


.  6:  Statistical  analysis  diagram  of  floating  point 
realization  method 
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This  result  is  better  than  the  output  signal-to-noise  ratio 
based  on  the  FFT  fixed  point  realization  and,  indeed,  this  i3 
because  the  operation  error  is  ignored. 

So  far  we  have  carried  out  a  detailed  analysis  of  the 
characteristics  of  DFT  of  DM  and  compared  the  result  with  that  of 
FFT.  In  the  course  of  our  analysis  all  the  inputs  were  assumed  to 
be  the  white  noise.  In  reality,  delta  modulation  can  not  be 
carried  using  this  type  of  signal .  Computer  simulation  shows  that 
if  the  input  were  the  sinusoidal  signal,  under  the  condition  that 
the  minimum  value  is  not  lost  when  the  summation  is  calculated; 
namely,  when  bq06,  br>-18,  the  increase  in  one  digit  of  bq  will 
result  in  the  increase  of  signal-to-noise  ratio  by  6dB.  This  is 
consistent  with  the  theoretical  analysis  that  when  the  decimal 
point  of  bq  in  increased  by  one  digit  the  signal-to-noise  ratio 
will  increase  by  6dB. 

V.  Hardware  Realization  of  Delta  Modulation  DFR 

In  the  previous  section,  when  the  round-off  effect  of 
floating  point  realization  was  analyzed,  a  constant  C  was  not 
used  to  prevent  overflow  of  the  coefficient.  In  the  following 
analysis,  the  floating  point  realization  is  analyzed  by  equation 
(11),  similar  to  the  fixed  point  realization  case.  The  discussion 
will  be  focused  firstly  on  the  hardware  realization  of  floating 
point,  then  the  determination  of  constant  C,  and  finally  the 
hardware  realization  of  fixed  point. 

As  shown  in  figure  7,  the  DM  series  (+1,-1)  is  first 
transformed  to  the  binary  coding  (0,1)  and  then  input  to  the 
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key: 


Fig.  7 


1  -  counter 

2  -  clock 

3  -  adder 


Hardware  realization  diagram  of  floating  point 
calculation 
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controlling  end  of  the  adder  for  addition/subtract  ion  operation. 
The  clock  is  used  as  the  counter  input  and  the  counting  is  used 
as  the  ROM  address  code.  Each  coefficient  is  stored  in  ROM 
sequentially.  Addition/subtraction  operation  is  executed  if  the 
binary  code  is  0/1.  After  N  addition/subtraction  operations,  the 
output  is  X(k) . 

Equation  (11)  can  be  re-written  as  follows.  Let 
g(n,k)«  © 

v(t)  -  CD 

■  - o 

According  to  the  cyclic  behavior  of  g(n,k),  the  above  equation 
can  be  expressed  as : 


X({)  “  X  dx(n)g(t>,  O  +  yi  Jx(n)g(ny  j) 
a«0  nmN/i 

Nil- 1 

—  5]  dx(n)g(n,  {)  -+•  (— l)*  s  ix  ("  +  O.  (22) 


•  »0 


The  above  equation  has  the  property  that  if  k=odd,  then  the 
N/2's  dx(n)  after  n>-N/2  have  opposite  signs  while  the  value  of 
g(n.k)  is  unchanged;  if  k-even,  then  the  N/2’s  dx(n)  after  n>-N/2 
have  their  signs  unchanged.  Combining  these  two  conditions,  the 
ROM  storage  space  can  be  reduced  by  half  and  only  N/2  is  needed 
for  maximum  counting.  When  n>-N/2,  the  counter  will  return  to  0 
and  ROM  will  be  re-read;  only  when  k  is  odd  a  pulse  will  be 
output  after  the  counter  is  full  to  control  the  input  binary  data 
to  transform  0  to  1  and  vice  veers. 

Similarly,  if  the  value  of  k  is  taken  to  be  a  factor  of  N, 
equation  (22)  can  be  further  decomposed  to  2k  terms  and  the  ROM 
storage  can  be  reduced  to  N/2k. 
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The  determination  of  C  of  g(n,k)  of  equation  (21)  can  be 


divided  into  the  real  and  imaginary  parts  since  g(n,k)  is  a 
complex  function.  Let 

«nt»(IRcU(-*,0)l.  -2-*-),  (23) 

min(|Rc(*(*,*))l,  |lm(f(»,O)|)-AC-2'O.*,"’2-,(0  (24) 

where  bc  is  the  length  of  the  valence  code  and  bm  is  the  length 

of  the  tail  number.  Solving  for  the  above  two  functions,  one 

would  obtain 


c 


-L  /i 
Ay  N 


(1 


-  2'*-). 


(25) 


The  value  of  bm  can  be  determined  if  the  signal-to-noise 
ratio  is  assigned.  C  can  then  be  calculated  using  the  above 
equation.  As  a  result.  bc  can  be  determined. 

In  the  above  analysis  the  word  length  of  g(n,k)  is 
determined.  Actually,  the  word  length  of  ROM  is  determined. 
Similarly,  the  word  length  of  the  adder  can  be  determined.  In 
order  to  ensure  that  the  minimum  quantity  of  g(n,k)  will  not  be 
lost  when  summation  is  calculated,  the  word  length  should  be 
determined  based  on  the  maximum  value  of  X(k)  and  the  minimum 
value  of  g(n,k).  Because  of  the  greater  difference  between  the 
maximum  value  of  X(k)  and  the  minimum  value  of  g(n,k),  the 
required  tail  number  is  longer. 

Figure  8  shows  the  hardware  realization.  Except  for  extra 
adders  and  digit-shifting  temporary  storers,  the  other  parts  are 
similar  to  the  floating  point  hardware  realization.  If  the 
summation  is  calculated  by  dividing  the  series  not  into  pairs  of 
numbers  but  into  groups  of  2s  numbers,  then  each  number  needs  to 
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key: 


Fig.  8 


1  -  counter 

2  -  digit-shifting  storer 

3  -  adder 

4  -  digit-shifting  storer 

5  -  adder 

:  Hardware  realization  diagram  of  fixed  point 
calculation 
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be  shifted  by  s  digits.  If  the  maximum  of  s  is  taken  as  log2N, 
only  one  adder  is  required.  However,  the  digits  of  the  adder  is 
also  increased  correspondingly. 


VI.  Computer  Simulation 

In  order  to  illustrate  the  feasibility  of  delta  modulated 
DFT,  we  have  chosen  two  different  sinusoidal  signals:  (1)  the 
input  is  the  three  sinusoidal  signals  with  different  frequencies, 
and  (2)  the  input  signal  is  the  band-limited  Gaussian  signal,  for 
simulation  on  VAX11/780  computers. 

For  the  first  kind  of  signals,  three  frequencies:  1kHz, 
0.5kHz,  0.25kHz  are  assumed  with  magnitude  of  5V,  IV.  and  3V, 
respectively. 

Based  on  the  Niquest  sampling  criterion,  l/T-2kHz,  if  the 
required  input  signal-to-noise  ratio  is  taken  as  40dB  then  K  has 
the  value  of  32  and  the  delta  modulation  sampling  frequency 
f s-l/dt~K/T«64kHz.  Therefore,  let 


*(/)  —  10$in2«/|(  ■+•  5  sin  2*fi»  +  3sin2ir f,i. 


(26) 


where  fi«0.5kHz.  f2“lJ<Hz,  f3~0. 25kHz.  In  order  to  avoid  slope 
overload,  let 


A  >  1 

A'  (27) 

then  ^>-1.055  and  A  is  taken  as  1.055.  Figure  9  shows  the 
Fourier  spectrum  based  on  the  delta  modulation  method.  It  can  be 
seen  from  figure  9  that,  because  of  the  sinusoidal  signal  of 
limited  length,  the  three  partitions  of  the  homogeneous  wave  of 
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fg  occurs  at  0.75kHz.  Figure  10  shows  the  effect  on  signal-to- 
noise  ratio  when  bq  and  br  are  with  different  word  lengths.  When, 
br>-18  and  bq-4,  5,  6,  7,  the  computer  simulated  result  is 
greater  than  the  theoretical  predicted  result  by  lOdB.  This  is 
because  in  the  theoretical  analysis,  the  input  is  assumed  to  be 
the  white  noise  signal  while  in  the  simulation  the  input  is  the 
sinusoidal  signal.  However,  we  can  also  learn  that  when  br>“18, 
each  increment  in  bq  will  result  in  the  increase  in  signal-to- 
noise  ratio  by  6  to  lOdB,  which  is  very  close  to  the 
theoretically  predicted  result. 

When  the  input  is  Gaussian  signal  the  optimal  step  size  is 
difficult  to  determine.  An  empirical  formulation  for  the 
determination  of  optimal  step  size  is  given  in  reference  [8] . 
Figure  11 (a?  is  the  Gaussian  signal  spectrum  of  limit  fm-lkHz 
based  on  delta  modulation  DFT  and  11(b)  shows  the  Gaussian  signal 
spectrum  of  limi  1kHz  based  on  ordinary  DFT. 

VII.  Concluding  Remark 

The  multiplier-free  DFT  based  on  DM  proposed  in  this  paper 
has  similar  characteristics  as  the  FFT.  Because  of  the  major 
parts  used  are  ROM  and  adder,  the  hardware  realization  is 
relatively  easy.  The  theoretical  analysis  and  computer  simulation 
provided  consistent  results.  Because  of  the  limitation  of  ROM 
access  speed,  it  can  only  be  used  in  processing  mainly  speech  and 
seismic  signals. 
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Fig.  11:  (a)  Limited  1kHz  Gaussian  signal  obtained  through 
delta  modulation  DFT  and, 

(b)  Limited  1kHz  Gaussian  signal  obtained  through 
ordinary  DFT 
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